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ABSTRACT 

Motivated by an interpretation of z ~ 3 objects by Jimenez and Haiman (2006), we examine 
processes that control the fraction of primordial {Z — 0) gas, and so primordial stars, in high-SFR 
Lyman break galaxies. A primordial fraction different from 1 or requires microscopic diffusion 
catalyzed by a velocity field with timescale comparable to the duration of star formation. The only 
process we found that satisfies this requirement for LBGs without fine-tuning is turbulence-enhanced 
mixing induced by exponential stretching and compressing of metal-rich ejecta. The time-dependence 
of the primordial fraction for this model is calculated. We show that conclusions for all the models 
discussed here are virtually independent of the IMF, including extremely top-heavy IMFs. 
Subject headings: ISM: evolution-galaxies: abundances-galaxies: evolution- galaxies: ISM-turbulence 



1. INTRODUCTION 

Galaxies begin their lives with entirely primordial {Z = 
0) gas. As they age, metal production and mixing can 
only reduce the primordial gas fraction. We have ex- 
plored the expected time dependence of the primordial 
fraction for various types of mixing and chemical evolu- 
tion processes using a kinetic equation for the evolution 
of the abundance distribution. The present Letter ad- 
dresses the narrower question of whether any models for 
mixing and chemical evolution predict this transforma- 
tion occurs at an accessible redshift. 

Jimenez and Haiman (2006) (hereafter JH) showed 
that several UV properties of a variety of objects, mostly 
Lyman break galaxies (LBGs), at redshift z ~ 3, can 
all be understood if these objects contain a substantial 
fraction, about 10 to 50 percent, of massive stars with 
essentially zero metallicity {Z ^ 10~^ Zq; we use Z 
and metallicity indiscriminantly here). These UV prop- 
erties cannot together be explained by a top-heavy IMF, 
and require Z = stars; a top-heavy IMF is not re- 
quired. Massive stars have short lifetimes, so their metal 
abundances reflect that of the concomitant gas. Thus, 
if JH are correct, a substantial fraction of the interstel- 
lar medium of these galaxies, with star formation ages a 
few hundred million years, has not been polluted by any 
products of nucleosynthesis.^ 

Motivated by JH and other suggestions for Z = stars 
(|Malhotra fc Rhoads 200^ : IShimasaku et al. 20061 ). and 
the fact that during some early period in the lives of 
all galaxies a transformation from primordial to non- 
primordial must occur, leaving spcctrophotometric sig- 
natures (JH; sec Schaerer 2003), we examined the viabil- 
ity of a number of models for this transformation. 

A change in the primordial fraction requires micro- 
scopic diffusion enhanced by a complex velocity field, 

^ The existence of Z = gas requires an IMF deficient in 
stars with M < 1 Mq. Otherwise the number of low- mass stars 
observable today would be large, contradicting observational limits 
on the star fraction with very small Z (see Oey 2003) by several 
orders of magnitude. We discuss the effects of different IMFs below, 
where we show that all our conclusions are independent of the IMF 
as long as it satisfies the requirement M > 1 Mq. 



such as instabilities in swept-up supernova shells or a 
turbulent interstellar medium (ISM), along with disper- 
sal of nucleosythesis products over many kpc, involving a 
scale range of ~ 10^. Existing hy drodynamic simulations 
in a cosmological conte xt (e.g., iGovernato et al. 2004 
IScannapieco et al. 20051 ) therefore cannot represent the 
transition of the primordial gas or mixing in general, 
since they adopt mixing rules that arbitrarily spread met- 
als among nearest neighbor cells or SPH particles. The 
only hydrodynamic simulations of true mixing of tracers 
in galaxies were concerned either with turbulent disper - 
sal of mean metallicity, not mixing (iKlessen fc Lin 2003(). 
or mixing of initial spatially periodic inhomogeneities by 
numerical diffusivity in a turbulent galaxy with no co n- 
tinuing source of metals (jde Avillez fc Mac Low 20021 ). 

Instead of simulations, the arguments given here are 
phenomenological, in order to clarify the essentials for 
each physical process. A detailed discussion is given 
elsewhere (Pan and Scalo 2006, hereafter PS) using a 
formal probability distribution evolution equation. Here 
we report that most processes we examined are either 
far too fast or slow to partially erase the primordial frac- 
tion (Sec. 2), except for mixing by turbulence-enhanced 
diffusivity based on exponential stretching of blobs of 
nucleosynthesis products, discussed in Sec. 3.^ 

2. ESTIMATES OF TIMESCALES FOR MIXING 
PROCESSES 

2.1. Star Formation Age 

The existence of a primordial gas fraction P{t) that 
is not nearly unity or zero, i.e., both P{t) and 1 — P{t) 
are significantly larger than zero, a condition we refer to 
as a "significant" or "intermediate" primordial fraction, 
requires a mixing or depletion process whose character- 
istic timescale is comparable to the star formation (SF) 
age, the time since SF began. If the mixing timescale 
is much smaller, P{t) will be nearly zero; if it is much 

A recent estimate of the effect of SN mixing (sec. 2.4 below) 
on the primordial fraction at high redshift (Tumlinson 2006) ap- 
parently used a shell mass much larger than found in analytical 
and numerical calculations of supernova remnant evolution (see 
[Thornton et al. 1998IJHanavama fc Tomisaka 20061) . 
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larger, P{t) will remain near unity. SF ages for LBGs 
and likely-related objects at somewhat different redshifts 
have been estimated by Papovich et al. (2001), Shapley 
et al. (2001), Erb et al. (2006) and others, using galaxy 
evolution models that assume an exponentially decreas- 
ing SFR that began at some time in the past. Although 
there is much variation, median SF ages are 3 x 10^ yr. 

These SF ages cannot be significantly different for a 
number of reasons. The lookback times of about 11.5 
Gyr for z =3 imply a strict upper limit to the SF age of 
2 Gyr, and a more likely upper limit of 1.5 Gyr (corre- 
spondi ng to z ^ 8). T he irregular morphologies of the 
LBGs (jGiavahsco 20021) suggest that LBGs are in the 
process of formation, accum ulating large fra gments of gas 
and stars through mergers (IConselice 20061 ). and proba- 
bly undergoing one of their first major bursts of SF; star- 
burst po£ulationsJiave_d^ statis- 
tics ()Kennicutt et al. 1987t iNiko lic et al. 2004'). model- 
ing of integrated light features (Marcillac et al. 2006), 
and theoretical arguments fsee [Leitherer 2001i) . that are 
similar to the estimates in LBGs. Finally, ages much 
greater than 3 x 10* y r would produce greater metallic- 
ity than observed (see iGiavahsco 200^ . 

We assume star formation has been a continuous func- 
tion of time. If instead the SFR preceding or during the 
present episode consists of bursts of shorter duration, 
most of our arguments remain unchanged if the SF age 
is replaced by the accumulated duration of SF. 

2.2. Depletion of the Primordial Fraction by Sources 

The JH result of 10-50% primordial at z ^ 3 seems 
surprising, but actually galaxies should remain almost 
completely primordial for billions of years in the absence 
of microscopic diffusivity. SN metal production slowly 
depletes primordial gas by transferring it from a, Z = 
delta function in the Z probability distribution (pdf) to 
another delta function at a much larger Z, the source 
metallicity Zs averaged over the IMF (~ 0.1 assuming 
the hot ejecta are well mixed). Intermediate values of Z 
cannot be reached without diffusivity. 

This suggests the simplest explanation for an inter- 
mediate primordial fraction in LBGs: 50 to 90% of the 
primordial gas passed through stars that became SNe. 
The timescale for this process is the source timescale, 
Tsrc = Mgas/BRsN, where B is the star formation 
rate and Rsn is the returned fraction from SNe aver- 
aged over the IMF. We take the total gas mass Mgas as 
5 X 10^^° Mq, extrapolated from gas masses estimated 
in the z ^ 2 UV-selected sample of Erb et al. (2006). 
From a number of studies, we adopt a median SFR of 100 
Mq/ji for L BGs at z ^ 3 a ssuming the IMF lower limit 
Ml = 0.1 (iPapovich e t al. 20 011: IShaplev et al 20011 : 
iGiavalisco 20021 : lErb et al. 200QnYan et al. 2006f) . For 
this Ml, we calculate Rsn ~ 0.1 for various IMFs using 
the ejected masses given in Woosley & Weaver (1995), 
Meynet & Maeder (2002) and Nomoto et al. (2006), find- 
ing little dependence on metallicity, including Z = 0, and 
20-30% variation between studies. Variations in the form 
of the IMF change the SFR by ~50%, with only a slight 
effect on Rsn- The source timescale is then Tsrc = 5 
Gyr, within a factor of a few considering the uncertainty 
in the SFR and Mgas, too large to affect the primordial 
fraction. 

We examined the IMF-dependence of Tsrc- There are 



two nonstandard IMFs that are especially relevant. 1. 
Intermediate primordial fractions require that the IMF 
lower limit M; ^ 1 Mq to avoid too many Z = stars 
observable today. Such a cutoff does not affect the source 
timescale: The empirical SFRs are based on integrated 
light from massive stars corrected for the rest of the IMF, 
and the resulting decrease of the SFR due to the cutoff is 
exactly compensated by the increase of the mass ejected 
by supernovae per unit mass of stars formed Rsn- This 
IMF-independence of Tsrc holds for any cutoff smaller 
than the lower mass limit for SNe, ~ 8 Mq. By the same 
argument, such a cutoff does not overproduce metallicity. 
2. A perennially popular IMF for Z = star formation 
consists of only very massive stars (VMS) due to the 
Jeans mass resulting from H2 cooling (Hutchins 1976; 
see Bromm and Larson 2004) altho ugh it has been ques- 
tioned on a number of grounds (|Silk fc Langer 2 006). 
Comparing Ha emission per unit SFR for a VMS IMF 
(50-500) Mq of Z = stars in Schaerer (2003) with 
the same quantity for a (0.1-100) Mq IMF in Kennicutt 
(1998), both for a Salpeter IMF (for illustration only), 
the SFR for a VMS IMF is 26 times smaller. Assuming 
only stars in the range 130-260 M p, explode as pair in- 
stabihty SNe (jWooslev et al. 2002D . we find Rsn = 0.27, 
so Tsrc = 50 Gyr, an order of magnitude larger than the 
normal IMF case. 

These results strengthen our conclusion that the forma- 
tion of massive stars is far too slow to deplete primordial 
gas significantly in the available time. 

2.3. Filling the Gap by Diffusion 

Without a process to spread metals into the "gap" be- 
tween the Zs peak and the primordial Z = peak, P{t) 
would remain near unity for billions of years. The only 
physical process that can fill this gap is microscopic diffu- 
sivity. However an estimate of the rate at which diffusiv- 
ity from random sources could pollute primordial gas in 
LBGs, using diffusion lengths for the cold neutral, warm 
neutral, and warm ionized ISM similar to those in Oey 
(2003), shows that the fraction of the mass of a galaxy 
mixed in time t is only ^ 10~^(i/0.5 Gyr)^/^, so diffu- 
sivity by itself cannot pollute more than a tiny fraction 
of the primordial gas over the estimated SF ages. 

To reduce the primordial fraction, a velocity field is 
required to catalyze diffusivity. However a velocity field 
cannot by itself affect the global metallicity distribution 
or primordial fraction, or mix at all: Displacement of 
fluid parcels of different Z by the velocity field conserves 
their volumes and thus volume fractions (or mass frac- 
tions for compressible flows), replacing one by another 
in space, having no effect on the Z distribution. This 
can be shown rigorously using a metallicity pdf equation 
(PS). A velocity field can only enhance mixing by spa- 
tially ramifying the Z field for diffusion to operate on 
small scales. Models that mix by sweeping of gas by SN 
or SB shells, cloud motions, differential rotation, or "tur- 
bulent diffusion" are unphysical without recognition that 
they are implicit models for microscopic diffusion. 

2.4. Several Catalyzing Velocity Fields 

Expanding SNRs and superbubbles (SBs) are the main 
agents of mixing in many sequential enrichment inho- 
mogeneous chemical evolution models (Reeves 1972; see 
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Tsujimoto et al. 1999; Oey 2000; Argast et al. 2000; 
Saleh et al. 2006) . Shells can mix, but only if 
instabilities allow diffusion to mix swept-up gas with 
new products of nucleosynthesis. Assuming this oc- 
curs, each SNR sweeps up and mixes a mass Msw 
2 X 10^ Mp, for Z = gas (jThornton et al. 19981 : 
iHanavama fc Tomisaka 200^ . The time to sweep up 
primordial ISM is Tsw = Mgas/{i^SNMsw), where the 
SN rate vsn = iB/{M^,). For an IMF with mass range 
(0.1, 100) Mq and indices (-0.4, -1.7) below and above 
1 Mq the number fraction e of stars that become SNe 
is 0.004 and the average stellar mass (M*) = 0.6, so 
vsN = 0.7/yr and = 3 x 10^ yr. Equivale ntly, the ac - 
cumulated volume fiUing factor NQ = t/Tsw (|Oev 20001 ) . 
so P = exp{—t/Tsw) = exp{—NQ) and it is impossible 
to preserve primordial gas in LBGs for longer than 1% 
of the observed SF ages ~ 3 x 10^ yr, unless the mixing 
efficiency is artificially tuned to 1%, in which case the 
model predicts too large a present-day scatter in metal- 
licity compared to observations. Superbubbles have a 
smaller frequency, but sweep up more gas, producing an 
almost identical result. Spatial clustering and infalling 
Z = gas also do not affect the conclusion. These points 
are discussed in details in a separate publication (PS). 
Using the same argument as for Tsrc in sec 2.2, an IMF 
cutoff at 1 Mq does not affect the SN rate, giving the 
same Tsw Unexpectedly, Tgw for a VMS IMF is also 
nearly unchanged, because of the cancellation of the de- 
crease of vsN (by a factor of 100) and the increase of 
Msw by a factor of ~ 50 due to the large explosion en- 
ergy of pair instability SNe (up to 10^^ erg, see Wooslcy 
et al. 2002). Therefore our conclusion that SN or SB 
sweep-up mixing is so fast that LBGs should have zero 
primordial gas is virtually independent of the IMF. The 
mixing would be 100 times faster using the mixed mass 
per event adopted by Tumlinson (2006). 

Unmixed pockets of metals in SBs could blast out of 
a galactic disk, later showering the disk with "droplets" 
of pu re metals, diffusively mixing once they land in the 
disk (iTenorio-Tagle 19961 ). Fine-tuning of the number 
of droplets, or equivalently the mixed mass, per SN, is 
required to give the desired timescale, but is unspecified 
by the model. 

Differential rotation could stretch the products of nu- 
cleosynthesis deposited in a ~ 100 pc SN blob, into long 
thin annuli until the scale of diffusivity is reached. The 
shear rate in LBGs, or whether they differentially rotate, 
is unknown. We used the rate ~ 10 km/sec/kpc in our 
Galaxy as an illustration and found that the timescale to 
reach the diffusive scale derived in sec. 3 below, is about 
8 Gyr, which is too slow. 

Another stretching process, turbulence, produces ex- 
ponential, rather than linear, stretching, with strain rates 
ten times larger than Galactic differential rotation. 



3. TURBULENCE-ENHANCED DIFFUSIVE MIXING 

Turbulence deforms large-scale fluid elements and the 
tracers they contain into filaments and sheets, bring- 
ing tracers closer together until scales are reached on 
which diffusivity can homogenize faster than the strain 
timescale of the fluid. The process is described by the 
general equation for the evolution of the metallicity field 



in an arbitrary velocity field u(x,<), 

^^^^+u(x,i)-VZ(x,i) = -V-{pKVZ{jc,t)) + S (1) 
ot p 

where k is the diffusivity and S denotes the sources. 

In our model a straining event on the scale of 
the sources Lg takes Z to a critical scale, Ldi/f, 
small enough for diffusivity to operate, at con- 
stant mean strain rate in a single step, by expo- 
nential stretchin g of line elements (|Batchelor 19521 : 
iVoth et al 2002^. This "short circuit" of the scalar 
cascade dVfllermau x et al. 20011) is support ed by exper- 
iments (Villcrmaux 2004; Voth ct al. 200 ^ and simula- 
tions ( Girimaii fc Pop e 1990: Goto fc Kida 2003[ ). and is 
similar to the scalar turbulence field theory of Shraiman 
and Siggia (2000). We assume that in supersonic tur- 
bulence compressions are analogous to stretching in the 
sense of bringing tracers to the critical scale Ldi / / . 

The scale L^iff below which the diffusivity term ex- 
ceeds the advection term in eq 1 is obtained by equating 
the two terms and replacing spatial derivatives by L^iff 
and u by the velocity at scale of L^iff. Assuming ui I 
scaling appropriate for exponential stretching, the result 
is Ldiff = [k/{U/Ls)]^^'^ , where U is the rms turbulent 
velocity on the scale of the sources. A residence-time 
average diffusivity, assuming the WNM contains more 
than 20% of the ISM mass, is k - 10^° cm^s'^ Then 
the average diffusivity scale in the ISM is Ldi/f ~ 0.06 

(Lioo/J7io)^^^ pc wher e Lipp is Ls/lOO pc, an d Uio is 
U in units of 10 km/s (jKulkarni fc Heiles 198% . We as- 
sume the same scale for z ~ 3 LBGs, noting that Ldiff 
only enters the mixing time (below) logarithmically. 

The mixing time follows from the assumed exponen- 
tial stretching, in which scales change as dl/dt = — 7Z, 
where 7 = U/Ls is the large scale strain rate. This 
gives a time to bring tracers from Ls to Ldi // Tmix = 
{Ls/U) ln{Ls/Ld^ff) = {Ls/2U) ln{ULs/K). The quan- 
tity ULs/k is the diffusivity analogue of the Reynolds 
number, called the Peclet number Pe. Numerically 
Pe 3 X 10'^ {U10L100/K20), giving r,„„ = {Ls/2U) 
ln{Pe) ^ 75 Myr. The result could be smaller if the 
velocity dispersion increas es with the SFR, as f or turbu- 
lence driven by SNe fe.g.. lDib fc Burkert 20051) . 

The primordial fraction is the cumulative pdf of the 

gas metaUicity f{Z' ,t)dZ' , where f{Z,t) is the dif- 
ferential pdf of metallicity, as the integration limit 
Z approaches zero. We can obtain the equation for 
P{t) heuristically, w ithout details of the integral closure 
panicka et al. 1979f) we used to derive the full pdf equa- 
tion corresponding to the advection-diffusion equation eq 
1 (PS), a turbulent mixing closure that gives the same 
timescale for exponential variance decay as found in sim- 
ulations by de Avillez and Mac Low (2002, eq. 17), and 
the same dependence of mixing time on initial size and 
SN rate (their Fig. 7), if velocity dispersion U scales as 
the square root of the SN rate (Dib and Burkert 2005). 

The primordial fraction decreases whenever ffuid ele- 
ments with primordial mass fraction P and gas that has 
been polluted by sources or previous mixing events, with 
mass fraction 1 — P, are stretched sufficiently to result in 
diffusive mixing. This interaction occurs with an average 
frequency t^]^- P is also gradually depleted by cycling 
through massive stars that inject metals when they ex- 
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SF duration (Gyr) 

Fig. 1. — Primordial fraction as a function of SF duration for 
combinations of tlie mixing timescale Tmix and the source timescale 
Tsrc (see sec 2.2). The empirical value for Tsrc is about 5 Gyr within 
a factor of a few for an IMF with M; = 0.1. Tsrc is independent of 
Ml if it is smaller than the mass limit for SNe. For a VMS IMF, 
Tsrc is much larger but does not affect the result significantly, as 
discussed in the text. The range suggested for 2 ~ 3 LBGs (JH) is 
indicated by the arrow. Also shown is an infall model with infall 
rate equal to the SFR. 

plode, on a timescale Tsrc(sec 2.2), which we assume is 
constant for the times of interest. The equation for the 
primordial fraction is then 

dP/dt = -P{1 - P)/t„„, - P/Tsrc (2) 

whose solution is, using the fact that Tmix ^ Tsrc, 

P = (1 + iT,mx/T'src)exp{t/T„rix)y^ (3) 

P decreases exponentially on timescale Tmix, but only 
after a delay time ~ Tmix ln{Tsrc/Tmix), which is 3 x 
10* yr for Tmix = 75 Myr and r^^c — 5 Gyr assuming an 
IMF with Ml <, 8 Mq. The delay time is the time for 
sources to provide enough non-primordial gas to make 
P depart from unity, but the dependence on Tsrc is only 
logarithmic, and so is nearly independent of IMF, even 
for the extreme case of a VMS IMF (see sec 2.2). 

The behavior of P as a function of time for different 
Tmix and Tsrc IS illustrated in Fig. 1 . The smallest mixing 
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timescale shown, 0.025 Gyr, corresponds to a larger ve- 
locity dispersion 30 km/s, not unreasonable for galax- 
ies with large SFRs. Our major result is that P{t) de- 
clines on a timescale similar to SF ages inferred from 
empirical modeling, without adjustment of parameters. 

The effect of Z = infall can be understood by adding 
to eq 2 a term (1 — P) /un where Tin = Mga^ /infall rate is 
the infall timescale. An example is shown in Fig. 1. Infall 
allows intermediate values of P(t) for a longer time, but 
only for infall timescales close to Tmix, implying a huge 
infall rate. Therefore it is unlikely that infall modifies 
the primordial fraction predicted by turbulent m ixing. 
Galactic winds with large rates (|Erb et al. 2006f) have 
no effect if the winds sample the full pdf of metallicity. 
If the galaxies have undergone previous episodes of SF 
with an accumulated duration as large as ^ 1 Gyr, even 
the turbulent model cannot explain the intermediate pri- 
mordial fractions claimed by JH. 

4. DISCUSSION 

Most mixing processes predict a primordial fraction 
that is either unity or zero at z ^ 3 because they mix on 
a timescale that is much larger or smaller than the empir- 
ical SF ages. P{t) should be zero in almost all galaxies if 
stellar explosions mix as efficiently as assumed in sequen- 
tial enrichment models (or much more efficiently, Tum- 
linson 2006). Our turbulence-enhanced diffusivity model 
naturally preserves primordial gas from rapid mixing for 
a few times the mixing time, which itself depends only 
weakly on parameters, in particular the assumed IMF or 
the averaged diffusivity, and gives an intermediate pri- 
mordial fraction in galaxies with SF ages ~ 1 — 3 x 10* yr. 
That this timescale happens to match the star formation 
ages of these galaxies is no coincidence if star f ormation 
is driven by turbulence (|Mac Low fc Klessen 2004) pow- 
ered by stellar explosions. Future systematic investiga- 
tions of the spectrophotometric signatures of primordial 
gas in galaxies could distinguish these possibilities. 
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